AFRL- VA- WP-TP-2006-3 1 0 

A  HYPERSONIC  VEHICLE  MODEL 
DEVELOPED  WITH  PISTON 
THEORY  (PREPRINT) 

Michael  W.  Oppenheimer 
David  B.  Doman 


JANUARY  2006 


Approved  for  public  release;  distribution  is  unlimited. 


STINFO  INTERIM  REPORT 

This  work  has  been  submitted  to  the  2006  AIAA  Atmospheric  Flight  Mechanics  Conference 
proceedings.  This  is  a  work  of  the  U.S.  Government  and  is  not  subject  to  copyright  protection  in 
the  United  States. 


AIR  VEHICLES  DIRECTORATE 

AIR  FORCE  RESEARCH  LABORATORY 

AIR  FORCE  MATERIEL  COMMAND 

WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7542 


NOTICE 


Using  Government  drawings,  specifications,  or  other  data  included  in  this  document  for  any 
purpose  other  than  Government  procurement  does  not  in  any  way  obligate  the  U.S.  Government. 
The  fact  that  the  Government  formulated  or  supplied  the  drawings,  specifications,  or  other  data 
does  not  license  the  holder  or  any  other  person  or  corporation;  or  convey  any  rights  or  pennission 
to  manufacture,  use,  or  sell  any  patented  invention  that  may  relate  to  them. 

This  report  was  cleared  for  public  release  by  the  Air  Force  Research  Laboratory  Wright  Site 
(AFRL/WS)  Public  Affairs  Office  (PAO)  and  is  releasable  to  the  National  Technical  Information 
Service  (NTIS).  It  will  be  available  to  the  general  public,  including  foreign  nationals. 

PAO  Case  Number:  AFRL/WS  06-0156,  18  Jan  2006. 


THIS  TECHNICAL  REPORT  IS  APPROVED  FOR  PUBLICATION. 


/s  / 

Michael  W.  Oppenheimer 
Electronics  Engineer 
Control  Design  and  Analysis  Branch 
Air  Force  Research  Laboratory 
Air  Vehicles  Directorate 


/s  / 

David  B.  Doman 
Acting  Chief 

Control  Design  and  Analysis  Branch 
Air  Force  Research  Laboratory 
Air  Vehicles  Directorate 


/s  / 

Brian  W.  Van  Vliet 
Chief 

Control  Sciences  Division 
Air  Force  Research  Laboratory 
Air  Vehicles  Directorate 


This  report  is  published  in  the  interest  of  scientific  and  technical  information  exchange  and  its 
publication  does  not  constitute  the  Government’s  approval  or  disapproval  of  its  ideas  or  findings. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  searching  existing  data 
sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis 
Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 


January  2006  Conference  Paper  Preprint  01/01/2005-01/09/2006 


4.  TITLE  AND  SUBTITLE 

A  HYPERSONIC  VEHICLE  MODEL  DEVELOPED  WITH  PISTON  THEORY 

5a.  CONTRACT  NUMBER 

In-house 

(PREPRINT) 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

N/A 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

Michael  W.  Oppenheimer  and  David  B.  Doman 

N/A 

5e.  TASK  NUMBER 

N/A 

5f.  WORK  UNIT  NUMBER 


_  |  N/A 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Control  Design  and  Analysis  Branch  (AFRL/VACA) 

Control  Sciences  Division 
Air  Vehicles  Directorate 

Air  Force  Research  Laboratory,  Air  Force  Materiel  Command 
Wright-Patterson  AFB,  OF!  45433-7542 _ 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Vehicles  Directorate 
Air  Force  Research  Laboratory 
Air  Force  Materiel  Command 
Wright-Patterson  Air  Force  Base,  OF!  45433-7542 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 

This  work  has  been  submitted  to  the  2006  AIAA  Atmospheric  Flight  Mechanics  Conference  proceedings.  This  is  a  work 
of  the  U.S.  Government  and  is  not  subject  to  copyright  protection  in  the  United  States. 

PAO  Case  Number:  AFRL/WS-06-0156  (cleared  January  18,  2006). 

14.  ABSTRACT 

For  high  Mach  number  flows,  M>  4,  piston  theory  has  been  used  to  calculate  the  pressures  on  the  surfaces  of  a  vehicle.  In 
a  two-dimensional  flow,  a  perpendicular  column  of  fluid  stays  intact  as  it  passes  over  a  solid  surface.  Thus,  the  pressure 
at  the  surface  can  be  calculated  assuming  the  surface  were  a  piston  moving  into  a  column  of  fluid.  In  this  work,  piston 
theory  is  used  to  calculate  the  rigid  body  forces,  moments,  and  stability  derivatives  of  a  hypothetical  hypersonic  vehicle. 
Only  longitudinal  motion  is  considered  in  this  case  and  lateral  motion  will  be  included  in  subsequent  work. 

15.  SUBJECT  TERMS 

Piston  Theory,  Hypersonic  Vehicles 

17.  LIMITATION  18.  NUMBER 
OF  ABSTRACT:  OF  PAGES 

SAR  12 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 
Michael  W.  Oppenheimer 
19b.  TELEPHONE  NUMBER  (Include  Area  Code) 

(937) 255-8490 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

Unclassified 

Unclassified 

Unclassified 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFRL-VA-WP -TP-2006-3 1 0 


10.  SPONSORING/MONITORING 
AGENCY  ACRONYM(S) 

AFRL/VACA 

11.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER(S) 

AFRL-VA-WP -TP -2006-3 1 0 


l 


A  Hypersonic  Vehicle  Model  Developed  With  Piston  Theory 


Midiael  W.  Oppenheimer  * 

David  B.  Doman  t 

Air  Force  Research  Laboratory,  WPAFB,  OH  45433-7531 


Abstract 

For  high  Mach  number  flows,  M  >  4,  piston  theory  has  been  used  to  calculate  the 
pressures  on  the  surfaces  of  a  vehicle.  In  a  two-dimensional  flow,  a  perpendicular  column  of 
fluid  stays  intact  as  it  passes  over  a  solid  surface.  Thus,  the  pressure  at  the  surface  can  be 
calculated  assuming  the  surface  were  a  piston  moving  into  a  column  of  fluid.  In  this  work, 
piston  theory  is  used  to  calculate  the  rigid  body  forces,  moments,  and  stability  derivatives 
of  a  hypothetical  hypersonic  vehicle.  Only  longitudinal  motion  is  considered  in  this  case  and 
lateral  motion  will  be  included  in  subsequent  work. 


HSV  Model 

Figure  1  shows  the  hypersonic  vehicle  considered  in  this  work.  The  vehicle  consists  of  4 
surfaces:  an  upper  surface  and  three  lower  surfaces.  All  pertinent  lengths  and  dimensions  are 
in  units  of  feet  and  degrees,  respectively.  The  goal  is  to  apply  piston  theory  to  determine  the 
pressure  distribution  on  the  surfaces  of  the  vehicle,  which,  in  turn,  can  yield  the  forces  and 
moments.  On  the  upper  surface,  the  surface  is  modelled  as  a  piston  moving  into  a  column  of 
fluid  that  has  the  properties  of  the  freestream.  On  the  lower  surface,  the  surface  is  modelled 
as  a  piston  moving  into  a  column  of  fluid  that  has  the  properties  of  the  fluid  behind  the 
oblique  shock.  The  pressure  on  the  face  of  a  piston  moving  into  a  column  of  perfect  gas  is 
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(1) 


where  the  subscript  ”oo”  refers  to  the  steady  flow  conditions  past  the  surface,  Vn  is  the 
velocity  of  the  surface  normal  to  the  steady  flow,  is  the  speed  of  sound,  and  P  is  the 
pressure.  Taking  the  binomial  expansion  of  Eq.  1  produces 

—  =  l+  27  7  —  1  Vn  =1  +  1 Yol 
Poo  7  1  2  ^oo  ^oo 
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Multiplying  through  by  P x  and  using  the  perfect  gas  law  (P  =  pRT)  and  the  definition  of 
the  speed  of  sound  (a2  =  7 RT)  yields  the  basic  result  from  first-order  linear  piston  theory 

P  Poo  ~t~  Poo®oobr).  (3) 

The  infinitesimal  force  due  to  the  pressure  is 

cJF  =  -PdAn  (4) 

where  dA  is  a  surface  element  and  n  is  the  outward  pointing  normal. 

To  compute  the  forces,  moments,  and  stability  derivatives,  consider  small  perturbations, 
from  a  steady  flight  condition  at  in  the  velocities  u,  v,  and  w  and  the  rates  p ,  g,  and  r. 
Consider  first  the  upper  surface.  The  velocity  of  a  point  on  the  upper  surface  due  to  these 
perturbations  is 

Vu  =  (Foo  cos(o;  +  3°)  +  u)  i  +  (VA,  sin(a  +  3°)  +  w)  k  +  u>  x  ru  (5) 

where  i,  k  are  unit  vectors  in  the  x  and  z  body  axes,  respectively,  is  the  angular  rate  vector 
and  a  is  the  angle  of  attack.  For  longitudinal  motion  only,  cj  =  qj  where  j  is  a  unit  vector  in 

the  y  body  axes  direction.  In  Eq.  5,  ru  is  the  position  vector  of  a  point  on  the  upper  surface 

given  by 

ru  =  ruxi  +  ruyj  +  ruzk  (6) 

For  the  lower  surface  defined  by  the  points  c  and  d  in  Figure  1,  use  the  velocity  of  the  flow 
after  the  oblique  shock  to  obtain 

Vied  —  (V2  cos  6.2°  +u)i  +  (~V2  sin  6.2°  +  w)  k  +  uj  x  ricd  (7) 

while  for  the  surface  defined  by  points  g  and  h 

Vigh  =  (V2  +  u)  1  +  (-F2  +  w)k  +  ux  rlgh  (8) 

where  r icd  and  r igh  are  position  vectors  of  a  point  on  the  lower  surface  given  by 


ricd  =  ncJ  +  ncdVj  +  rlcdZk 

(9) 

righ  =  righxi  +  rlghyj  +  rlghZk 

(10) 

An  expansion  fan  occurs  on  the  final  lower  surface  and  will  be  considered  shortly.  The  normal 
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vectors  for  the  upper  and  lower  surfaces  are 


^■u  Tlux^  “I-  Tluyj  “I-  'ft'uzk 

n/cd  =  nicJ  +  nicdyj  +  mcdZk  (11) 

nigh  =  nighxi  +  nlghVj  +  nXghZk 

Performing  the  cross  products  required  by  Eqs.  5,  7,  and  8  gives 

i  j  k 

cj  xru=  0  q  0  =  qruzi  -  qruxk  (12) 

r/y*  /p 

_  ux  1  uy  1  uz 

i  j  k 

U  x  rlcd  0  q  0  qrlcdZi  -  qrXcdXk  (13) 

.  TlcdX  rlcdV  rlcdZ  . 

i  j  k 

u  x  rlgh  0  q  0  qnghZi  -  qrXghXk  (14) 

.  rlghX  T-lghV  rlghZ  . 

To  determine  the  velocity  of  interest,  i.e. ,  the  velocity  of  the  surface  normal  to  the  steady 
flow,  take  the  dot  product  of  the  velocity  with  the  appropriate  normal  vector  and  use  the 
result  in  Eq.  3  to  obtain 

Pu  Poo  T  PoO®OC  (Vu  •  nu) 

Picd  =  P2  +  p2a2  (Vicd  •  nlcd)  (15) 

Pigh  =  P2  +  P2&2  (Vigh  ■  n igh) 

Substituting  the  results  of  Eq.  15  into  Eq.  4  gives 

dF  u  {  P x)  Poo&oo  ( dAuiiu 

dFicd  =  {-P2  -  P2a2  (Vicd  •  n,J}  dAicdnicd  (16) 

dFigh  =  {~p2  ~  P2CI2  (Vlgh  ■  nlgh)}  dAlghnlgh 

Using  (Eqs.  5,  7,  and  8)  and  the  appropriate  normal  vector  (Eq.  11),  Eq.  16  becomes 


dFu  =  (-Poo  -  Poofloo  {(Vo  cos  (3°  +  a)  +  u  +  qruz )  nux  +  (Uoo  sin  (3°  +  a)  +  w  -  qrux)  nuz })  dAunu 
dF icd  =  (~P2  -  p2a2  {(V2  cos  6.2°  +  u  +  qrtcdZ)  nXcdX  +  (V2  sin  6.2°  +  w-  qrXcdX)  nlcdZ })  dAlcdnXcd 

dFigh  =  (- p2  ~  P2CI2  {  (V2  +  u  +  qrlghZ )  n,cdX  +  (V2  +  w  -  qrlghX )  nXghZ})  dAlghnlgh 

(17) 
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The  upper  and  lower  surface  elements,  cL4jn,:  can  be  written  as 


dAunu  =  [nuxi  +  nuyj  +  nuzkj  dAu 
dAicdnicd  =  (nlcdJ  +  mcdyj  +  nicdZk^j  dAtcd  (18) 

dAiBhniBh  =  +  niBhvJ  +  niahzk)  dAlgh 

In  order  to  evaluate  the  forces,  the  required  integrations  must  be  performed.  In  order  to 
perform  this  task,  the  position  vectors  and  normal  surface  vectors  must  be  determined.  For 
the  upper  surface,  the  position  vector  can  be  computed  by  evaluating  the  geometry  in  Fig.  1 
to  yield 

r u  =  xi  +  tan  3°  (x  —  55)  k  .  . 

-45  <  x  <  55 

The  upper  surface  normal  vector  is 


hy  —  sin  ?Ai  +  cos  3  °k 


(20) 


For  the  lower  surface,  these  position  vectors  are 


rtcd  =  xi  +  {—  tan 6.2°  (x  —  8)  +  b}k 
8<x<55 


(21) 


I';  ,..  =  xi  +  (b  +  hi)  k 
-12  <x  <S 


(22) 


while  the  normal  vectors  become 


hicd  =  —  sin  83. 8°  i  +  cos  83. 8°  k 

« I,,,  = k 


(23) 


To  compute  the  aerodynamic  properties  of  the  vehicle,  the  quantities  in  Eq.  17  are  integrated 
over  the  surface  of  the  vehicle.  Using  Eqs.  17-  19,  force  on  the  upper  surface  is 


0  /-55 


Fu  = 


-l  J  - 45 


(—Poo  +  PooOooA)  (  —  sin  3 °i  +  cos  3 °k  )  dxdz 


(24) 


where 


A  =  {(V^o  (3°  +  a)  +  u  +  qx)  tan  3°  (x  —  55)  +  (I4c  sin  (3°  +  a)  +  w  —  qx)  cos  3°}  (25) 
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For  the  lower  surfaces,  the  forces  become 


0  r  55 


Ft,  = 


b  J  8 


(P2  +  p2a2B)  ( —  sin  83.8°f  +  cos  83. 8° k )  dxdz 


(26) 


where 


B  =  {V2  cos  6.2°  +  u  +  q  (—  tan  6.2°  {x  —  8)  +  b )  (—  sin  83.8°)  +  (V2  sin  6.2°  +  w 


and 


fibh=  f  {~P2  +  P2a2(w 
J- 12 


Performing  the  integrations  produces 


qx )}  kdx 


qx )  cos  83.8°} 

(27) 

(28) 
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Fu  =  100/Pqo  —  5000  ( V'oo  cos  (3°  +  a)  +  u)  la^p^  tan  3° 
+1000,00^00  (Voo  sin  (3°  +  a)  +  w)  l  cos  3° 

<W>oo^tan3°  —  500aoopoolqtan3°q  (—  sin  (3°  +  a)  i  +  cos  (3°  +  a)  k 


(29) 


Figh  —  20  (— P2  +  a2p2  +  2g})  k  (30) 

Ficd  =  47 P2b  —  47a2p2bV2  cos  6.2°  +  47a2p2bV2  sin  6.2°  cos  83.8° 

2209 

—A7a2p2b2V2  sin  6.2°  sin  83.8°  H - — — a2p2b tan.  6. 2°V2  sin  6.2°  sin  83.8°  (31) 

2961  /  - 

H — — — a2p2bq  cos  83.8°  sin  83.8 °i  +  cos  83.8°kJ 

Future  Work 

At  this  point,  the  forces  on  the  vehicle  have  been  determined.  These  will  be  used  to 
determine  the  pitching  moment  and  similar  analysis  will  be  used  to  evaluate  the  stability 
derivatives.  Also,  the  flexible  body  effects  on  the  vehicle  will  be  taken  into  account. 


Conclusions 

In  this  work,  piston  theory  is  used  to  develop  a  model  for  the  longitudinal  dynamics  of  a 
hypersonic  vehicle.  In  particular,  velocities  of  flow  normal  to  the  surface  of  the  vehicle  are 
used  in  a  first  order  piston  theory  framework  to  determine  the  pressures  on  the  surfaces  of 
the  vehicle.  The  pressures  are  then  integrated  over  the  body  to  determine  the  forces  acting 
on  the  vehicle.  Piston  theory  is  useful  here  because  it  allows  the  inclusion  of  the  unsteady 
aerodynamic  effects,  which  are  not  captured  using  other  techniques. 
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